Transition temperature of a dilute homogeneous imperfect Bose gas 



Peter Arnold 

Department of Physics, University of Virginia, P.O. Box 400714, Charlottesville, VA 22904-4"^^ 



o 
o 

a\ 

m 
> 

00 
(N 
(N 
m 
o 



o 



Department of Physics, 



Guy Moore 

University of Washington, Seattle, 
(February 1, 2008) 



Washington 98195-1560 



The leading-order eflect of interactions on a homogeneous Bose gas is theoretically predicted to shift 
the critical temperature by an amount ATc ~ #asc«-^''^?o from the ideal gas result To, where Osc is 
the scattering length and n is the density. There have been several different theoretical estimates 
for the numerical coefflcient We claim to settle the issue by measuring the numerical coefficient 
in a lattice simulation of 0(2) cj)"^ field theory in three dimensions — an effective theory which, as 
observed previously in the literature, can be systematically matched to the dilute Bose gas problem 
to reproduce non-universal quantities such as the critical temperature. We find ^ = f .32 ± 0.02. 



The phase transition temperature Tq of an ideal three- 
dimensional Bose-Einstein gas at fixed density is some- 
thing that every physicist learns to calculate in gradu- 
ate school, if not before. It is amusing that the first 
correction to that result, from arbitrarily weak interac- 
tions, is sufficiently challenging that there has not yet 
been theoretical agreement on its magnitude. It is under- 
stood that, in the weak interaction limit, the correction 
ATc = Tc — Tq behaves parametrically as 
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where c is a numerical coefficient. (A clean argument may 
be found in Ref. Q.) Here Osc is the scattering length, 
and the weak interaction (or dilute) limit is agcTL^^^ <C 1: 
that is, flsc is small compared to the typical separation 
between particles. (We will assume that the interaction 
is repulsive.) However, there has been little agreement 
in estimates of the coefficient c, a variety of which 
are shown in Fig. 0. Some of these estimates are adver- 
tised as rough, but others are not. The difficulty arises 
because the phase transition is second order, and per- 
turbation theory typically breaks down at second-order 
phase transitions: the physics that determines ATc is 
non-perturbative. As we shall briefly review, the prob- 
lem of finding ATc iu the weak interaction limit can be 
related to solving static three-dimensional 0(2) scalar 
field theory |]. In this paper, we present results from 
using standard, numerical, lattice Monte Carlo methods 
to solve that theory. In principle this provides an ex- 
act method for computing c to any desired precision; in 
practice, one is limited by computer time and memory. 
Working on desktop computers, we find c = 1.32 ± 0.02, 
which is the grey bar in Fig. ^. 

It is long distance physics that determines ATc in the 
weak interaction limit. It is well known |lO| ] that, at dis- 
tance scales large compared to the scattering length Ogc, 
an appropriate effective theory for a dilute Bose gas is the 
second- quantized Schrodinger equation, together with a 
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FIG. 1. Estimates from the literature of the constant c in 
ATc/Tq — > can^^^. The grey bar is the result of this paper. 



chemical potential /i that couples to particle number den- 
sity 4'*4': ^-nd a IV'I^ contact interaction that reproduces 
low-energy scattering. The corresponding Lagrangian is 



ir^?. (2) 



Corrections to this effective theory (due, for instance, to 
energy dependence of the cross-section or 3-body inter- 
actions) may be ignored for the propose of computing 
the leading-order result for ATc. To study (|^) at fi- 
nite temperature, apply the imaginary time formalism, 
so that t becomes ir and imaginary time r is periodic 
with period h(3 = fi/k-QT. The field ij] can then be de- 
composed into frequency modes with Matsubara frequen- 
cies LOn — 2Tm/hp. For distances large compared to the 
thermal wavelength A = ?i\/27r/3/m, and sufficiently near 
the transition so that \fj,\ <C T, the non-zero Matsubara 
frequencies decouple from the dynamics, leaving behind 
an effective theory of only the zero- frequency modes tAoi 
with the action S = dr J <Px C becoming [|lj 
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up to corrections that again do not affect the leading- 
order result for ATc. This action can be thought of as the 
PH of a classical three-dimensional field theory. Finally, 
it is convenient to rewrite ■00 = +i4>2)V2n/X so that 
the effective action becomes a conventionally normalized 
0(2) field theory: 
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where (p is understood to be a 2-componcnt real vector 
{4'ii4>2) and u — 967r^asc/A^. 

As noted by Baym et al. , it is technically somewhat 
more convenient, in this formalism, to calculate the shift 
AriciT) in the critical density at fixed temperature in- 
stead of the shift ATc(n) in the critical temperature at 
fixed density. The two are trivially related at first order 
in ATc by AT^/T = -|Anc/?i, where the factor of | 
arises from the ideal gas relation Tq (x n^/^. In the field 
theory, n is given by {'ip*^), which is proportional to ((/>^). 
Putting everything together 
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is the difference between the effective theory value of 
{(fP'), at the critical point, for the cases of (i) u small 
and (ii) the ideal gas u—Q. Unlike {(fP')c, the difference 
A(0^)c is an infrared quantity, independent of how the 
effective theory (0) is regularized in the ultraviolet (UV) . 

Finding [(</'^)cJ^ in the effective theory corresponds 
to fixing M, varying r to reach the critical point, and 
then measuring (0^). The only parameter of this prob- 
lem is u. By dimensional analysis, physics is therefore 
non-perturbative at the infrared length scale and, 
again by dimensional analysis, A(0^)c is proportional 
to u. Putting this together with the ideal gas formula 
To — (27r?i'^/fcBm)[n/C(|)]^/'^, one may summarize the 
relationship between the weak interaction limit for ATc 
and the 0(2) effective theory as ATc/Tq c a^c n^^^ with 
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To compute A(0^)c/w, we put the 0(2) theory (H) on a 
lattice. For example, the most straightforward discretiza- 
tion would use the action 
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on a simple cubic lattice, where (xy) represents all 
nearest-neighbor pairs and aiat is the lattice spacing (un- 
related to Osc)- The dimensionless coupling of the lattice 
theory is waiat, and the continuum limit is uaiat — > 0. 

Our simulations use an improved action to reduce lat- 
tice spacing errors. (The subscript U on S* in Eq. (||) 
stands for "unimproved.") Details concerning the action 
and our simulations are given in Ref. Our simu- 

lations use a combination of heat bath and multi-grid 
updates jl2|. At finite volume, we use the method of 
Binder cumulants |Tl|jr^ ] to determine a nominal critical 
value for r. 

Because three-dimensional scalar theory requires UV 
regularization of its cjP' interactions, the riat in the lattice 
action (B) is not simply the r^s in the effective theory 
action (4), which in turn is not simply related to the 
chemical potential /i in the original action (^. However, 
for our purposes here, we are only interested in adjusting 
riat to find the critical point, for a given u, and measuring 
A((/)^) there. The actual relation between riat, '"eff, and /i 
is unnecessary [|4|. The 0^ interactions, in contrast, do 
not require UV regularization: in the limit uaiat — > 0, the 
coupling u of the lattice theory may be identified with the 
continuum coupfing u = 967r^asc/A^ introduced earfier. 

For a given u, we compute (0^)c by Monte Carlo nu- 
merical simulations of the lattice theory. The w=0 piece 
of the difference A(0^)c can be easily computed without 
simulations: 
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where G[^j(p) is the free lattice propagator, and the mo- 
mentum p is integrated over the Brillouin zone (BZ). 
Such integrals are reviewed in Ref. For the unim- 

proved lattice theory (0), for example |16|, 
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where Sy = 3.175 911 535 625 • • •. 

There are two limits that must be taken of lattice 
Monte Carlo data: the contimmm limit uajat— >0 and 
the infinite volume limit Lu—>oo. Fig. || shows the de- 
pendence of our data on system size (Lu) at uaiat = 6. 
Fig. ^ shows the dependence on lattice spacing (woiat) at 
Lu = 576. From these two figures, it is reasonably clear 
that our raw data includes reasonably large volumes and 
reasonably small lattice spacings. We will discuss ex- 
trapolations of the infinite- volume continuum limit using 
finite volume scaling, the known critical exponents of this 
model, and an analysis of finite lattice-spacing errors. 

First, however, we wish to show that one can obtain 
a quick estimate from the raw data without relying on 
anything fancy. Table | shows values associated with 
a selected subset of the largest Lu and smallest uaiat 
data points from the figures. We take the (Lu,Maiat) = 
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FIG. 2. 0(itaiat)-corrected results for A{cj>^)c vs. system 
size at uaiat = 6. A numerical fit to large L scaling behavior is 
shown, which fits the last 4 points with confidence level 61%. 
Also shown, for comparison, is a small Lu expansion |l^ of 
the exact continuum result in finite volume (at next-to-leading 
order in Lu). 
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FIG. 3. Resuhs for A(<^2>c vs. tiOiat at Lu — 576 that in- 
corporate perturbatively calculated O(iiaiat) corrections. The 
line is a fit of j4 + -B(Maiat)^ to all but the rightmost data point 
and has confidence level 14%. 
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(576, 6) result as a starting point for our estimate. The 
finite-volume correction is at least as big as the differ- 
ence with the value at (1152,6) [which corresponds to 
doubling Lu] but is unlikely to be double this difference. 
This difference is roughly —0.00015 (ignoring the small 
statistical errors) , so we might estimate the finite- volume 
correction to the Lu — 576 value to be somewhere be- 
tween —0.00015 and —0.00030. From a similar compari- 
son of the (576,6) and (576,3) data, we might estimate the 
finite lattice spacing correction to (576,6) to be between 
roughly -(-0.00005 and -1-0.00011. Adding our corrections 
and the original (576, 6) data point, the final rough esti- 
mate of the continuum infinite-volume value would be 



-0.00119 ± 0.00011, 
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which, by (|^) would translate to c = 1.31 ± 0.12. 

We now summarize a more careful analysis of correc- 
tions, detailed in ref. [Q. Our strategy is to start again 
with Lu = 576 data, extrapolate the continuum limit, 
and estimate the finite volume correction. To improve 
the approach to the continuum limit, we have analytically 
calculated the O(uaiat) corrections to A(0^)c and the re- 
lation between lattice and continuum values of u, using 
lattice perturbation theory ||ll[. Fig. ^ shows Lu — 144 
data which clearly demonstrates our control of lattice 
spacing errors. The uncorrected data clearly has a linear 
dependence on uaiat- But the corrected data fits, to high 
confidence level, the assertion that the remaining error 
scales as (uaiat)^. Based on a similar fit to the data of 
Fig. 1^, we estimate the uaiat result at Lu — 576 as 
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Finite scaling arguments predict that the large Lu cor- 
rections to A((tP)c should scale as [0 
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if the method of Binder cumulants is used to determine 
the transition point in finite volume. Here a ~ —0.01 
and 1/ = (2 — Q!)/3 ~ 0.67 are the specific heat and cor- 
relation length critical exponents of the 0(2) model [pT| . 
Further discussions of fits, and an analysis of corrections 
to scaling, may be found in ref. ||ll[. On the basis of these 
fits, we estimate the finite size correction at Lu = 576 to 
be 0.000241 ± 0.000007 . Putting this together with the 
Lu = 576 continuum extrapolation (|l^), we obtain the 
final infinite- volume continuum result 
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-0.001198 ±0.000017. 
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Using Eq. (|^) for the weak- interaction limit of AT^, the 
result for A((j)'^)r translates to 



TABLE I. Selected data for A{<fP')^/u. 
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FIG. 4. The squares show A((;6^)c vs. uaiat at Lu = 144 
incorporating O(uaiat) corrections. The hne through them is 
a fit of the first 6 points to A + B(ua\^tf' and has confidence 
level 94%. The diamonds represent the corresponding uncor- 
rected data, with a straight line fit to the first 3 points to 
guide the eye. 



It is interesting to compare our numerical results with 
results from the large N expansion, depicted by triangles 
in Fig. 0. In large N ^ one generalizes the 0(2) effective 
field theory to an 0{N) theory of N real fields, calculates 
results in powers of 1/N, and hopes the expansion will be 
useful for the case of interest, N = 2. The leading-order 
(LO) result c ~ 2.33 [§ is off by 77%, but the next-to- 
leading-order (NLO) result c ~ 1.71 ||] moves in the right 
direction and is off by only 30%. This is surprisingly good 
for an expansion that treats = 2 as large. 

We should comment on the discrepancy of our results 
with previous numerical calculations in the literature 
1^,^, shown in Fig. ^ which used a radically different 
starting point. Rather than using field theory methods 
and the grand canonical ensemble, they start with the 
path integral for a large, fixed number of particles in a 
box. The recent work of Holzmann and Krauth , how- 
ever, makes a flawed assumption at the very beginning: 
they assume that the integrand of the path integral can 
be expanded perturbatively in the interaction, and keep 
only the leading term. This is wrong because interac- 
tions cannot, generically, be treated perturbatively at a 
second-order phase transition. 

We believe one likely problem with the older simu- 
lations of Griiter et al. is inadequate system size [ pT] . 
Reppy et al. have reported c = 5.1 ± 0.9 from exper- 
imental data on He-Vycor systems, but cautions about 
the data's interpretation may be found in Ref. ||^. 

As we completed this work, another paper appeared 
p9| which uses techniques very similar to ours and ob- 



tains the statistically compatible result c = 1.29 ± 0.05. 
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